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Abstract 

Dynamical analysis of spherically symmetric collapsing star sur¬ 
rounding in locally anisotropic environment with expansion-free con¬ 
dition is presented in f(R,T) gravity, where R corresponds to Ricci 
scalar and T stands for the trace of energy momentum tensor. The 
modified field equations and evolution equations are reconstructed in 
the framework of f(R,T ) gravity. In order to acquire the collapse 
equation we implement the perturbation on all matter variables and 
dark source components comprising the viable f(R,T ) model. The 
instability range is described in Newtonian and post-Newtonian ap¬ 
proximation. It is observed that the unequal stresses and density 
profile defines instability rage rather then the adiabatic index. How¬ 
ever, the physical quantities are constrained to maintain positivity of 
energy density and stable stellar configuration. 
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1 Introduction 


The astrophysics and astronomical theories are invigorated largely by the 
gravitational collapse and instability range explorations of self gravitating 
objects. Celestial objects tend to collapse when they exhaust all their nuclear 
fuel, gravity takes over as the inward governing force. The gravitating bodies 
undergoing collapse face contraction to a point that results in high energy 
dissipation in the form of heat flux or radiation transport [jT]. The end state of 
stellar collapse has been studied extensively, a continual evolution of compact 
object might end up as a naked singularity or as black hole depending upon 
the size of collapsing star and also on the background that plays important 
role in pressure to gravity imbalances 

The gravitating objects are interesting only when they are stable against 
fluctuations, supermassive stars tends to be more unstable in comparison to 
the less massive stars jo j. Instability problem in star’s evolution is of fun¬ 
damental importance, Chandrasekhar [5j presented the primary explorations 
on dynamical instability of spherical stars. He identified instability range 
of star having mass M and radius r by a factor T pertaining the inequal¬ 
ity T> § + n^-. Adiabatic index measures compressibility of the fluid i.e., 
variation of pressure with a given change in the density. The analysis of 
expanding and collapsing regions in gravitational collapse was presented by 
Sharif and Abbas [7]. 

Herrera and his collaborators 0-im presented the dynamical analysis 
associated with isotropy, local anisotropy, shear, radiation and dissipation 
with the help of T, it was established that minor alterations from isotropic 
profile or slight change in shearing effects bring drastic changes in range of 
instability. However, instability range of stars with zero expansion does not 
depend on stiffness of fluid, rather on other physical parameters [12]-pH], such 
as mass distribution, energy density profile, radial and tangential pressure. 
The impact of local anisotropy on plane expansion-free gravitational collapse 
is studied in [15] . 

General Relativity (GR) facilitates in providing field equations that leads 
to the dynamics of universe in accordance with its material ingredients. The 
predictions of GR are suitable for small distances, however, there are some 
limitations of GR in description of late time universe. Modified gravity theo¬ 
ries have been widely used to incorporate dark energy components of universe 
by inducing alterations in Einstein Hilbert (EH) action. Due to modifications 
in laws of gravity at long distances, dark source terms of modified gravity 
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leaves phenomenal observational signatures such as cosmic microwave back¬ 
ground, weak lensing and galaxy clustering [l]3j-[[2Ij]. Many people investi¬ 
gated the dynamics of collapse and instability range in modified theories of 
gravity, Cembranos et al. [22] studied the collapse of self-gravitating dust 
particles. Sharif and Rani |23j established the instability range of locally 
anisotropic non-dissipative evolution in f(T) theory. 

Among modified gravity theories, f(R) exhibits the most elementary 
modifications to EH action by adopting a general function f(R) of Ricci 
scalar. Ghosh and Maharaj [23] indicated that null dust non-static collapse 
in f(R) around de-Sitter higher dimensional background leads to naked sin¬ 
gularity. Combined effect of electromagnetic field and viable f(R) model has 
been investigated in |25| . concluding that inclusion of Maxwell source tends 
to enhance the stability range. Borisov et al. |26] investigated the spherically 
symmetric collapse of f(R) models with non-linear coupling scalar by execu¬ 
tion of one-dimensional numerical simulations. The dynamical instability of 
extremal Schwarzschild de-Sitter background framed in f(R) is investigated 
in [27] . 

Another modification of GR and generalization of f(R) was presented in 
2011 by Harko et al. [28] termed as f(R, T) gravity theory constituting the 
matter and geometry coupling, EH action is modified in a way that gravita¬ 
tional Lagrangian includes higher order curvature terms alongwith the trace 
of energy momentum tensor T. Shabani and Farhoudi |29| explained the 
weak field limit by applying dynamical system approach and analyzed the 
cosmological implications of /(/?,. T ) models with a variety of cosmological 
parameters such as Hubble parameter, its inverse, snap parameters, weight 
function, deceleration, jerk and equation of state parameter. Ayuso et al. 
[30] worked on consistency criterion for non-minimally coupled class of mod¬ 
ified theories of gravity. Sharif and Zubair [31] - [333 ascertained the laws 
of thermodynamics, energy conditions and analyzed the anisotropic universe 
models in f(R,T) framework. 

Chakraborty [35] explored various aspects of homogeneous and isotropic 
cosmological models in f(R, T) and formulate the energy conditions for per¬ 
fect fluid. Dynamics of scalar perturbations in f(R,T ) is explored in [35J. 
Jamil et al. [37] [38] reconstruct some cosmological models and studied laws 
of thermodynamics in f(R,T). In a recent paper [39], dynamical instability 
of isotropic collapsing fluid in the context of f(R, T) is considered. We have 
also discussed the stability analysis of spherically symmetric collapsing star 
surrounding in locally anisotropic environment in f(R,T) gravity [3D]- Fur- 
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thermore, conditions on physical quantities are constructed for Newtonian 
and post-Newtonian eras to address instability problem. 

Herein, we intend to develop the instability range of f(R, T) model under 
anisotropic background constraining to zero expansion. The expansion-free 
condition necessarily implies the appearance of cavity within fluid distribu¬ 
tion that might help in modeling of voids at cosmological scales. Also, such 
distributions must bear energy density inhomogeneities that are incorporated 
here by inducing non-constant energy density and pressure anisotropy. The 
dynamical analysis of various fluid distributions with expansion-free condi¬ 
tion has been studied in f(R) [iTj-ji3], however, such situations have not 
been covered yet in f(R,T). Recently, Noureen and Zubair [Hj discussed 
the implications of extended Starobinsky Model on dynamical instability of 
axially symmetric gravitating body. 

To develop collapse equation in f(R, T), we construct corresponding held 
equations constituting expansion-free fluid. The action in f(R, T) is as in 

+ £«]. ( 11 ) 

where C( m ) is matter Lagrangian and g denotes the metric tensor. The 
Lagrangian £( m ) can assume various choices, each choice corresponds to a 
set of held equations for some special form of fluid. Here, we have chosen 
£( m ) = p, 8 ttG = 1 and upon variation of above action with metric g uv the 
held equations are formed as 

(/t + 1)2~^ ~~ PQuvfr H - g U v 

+ (V U V„ — g uv O)fii] , (1-2) 

where denotes the energy momentum tensor for usual matter. 

The matter Lagrangian is configured in a way that it depends only on the 
components of metric tensor [45j. In order to present the dynamical anal¬ 
ysis we implement the linear perturbation on collapse equation, assuming 
that initially all physical quantities are in static equilibrium. The paper is 
arranged as: Einstein’s held equations and dynamical equations for f(R, T) 
are constructed in section 2 that leads to the collapse equation. In section 

3 perturbation scheme is implemented to the dynamical equations. Section 

4 covers the discussion of expansion-free condition and the components af¬ 
fecting the stability of gravitating objects, extracted from perturbed Bianchi 
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identities along with corrections to Newtonian and post-Newtonian eras and 
GR solution. Section 5 comprises the summary followed by an appendix. 


2 Dynamical Equations in f(R, T) 

We choose a three dimensional external spherical boundary surface £ that 
pertains two regions of spacetime termed as interior and exterior regions. 
The line element for region inside the boundary £ is of the form 

ds 2 __ = W 2 (t, r)dt 2 — X 2 (t, r)dr 2 — Y 2 (t, r){d6 2 + sin 2 ddcj?). (2.3) 

The domain beyond (lying outside) £ is exterior region with following line 
element [43j 

ds 2 + = ^1-—^ du 2 + 2 drdu — r 2 {d0 2 + sin 2 9d(p 2 ), (2.4) 

where v is the corresponding retarded time and M is total mass. To ar¬ 
rive at onset of field equations given in Eq. (11.21) . we choose describing 
anisotropic fluid distribution of usual matter, given as 

= (p + p±)V u V v - p±g uv + (p r - P±)XuXv , (2.5) 

where p denotes energy density, V u is four-velocity of the fluid, Xu corresponds 
to radial four vector, p r and p± represent the radial and tangential pressures, 
respectively. The physical quantities appear in energy momentum tensor are 
in accordance with the following identities 

V V = W~ 1 8Z, yU Vu = 1, X U = W-X, x u Xu = -l. (2.6) 

The expansion scalar 0 defines rate of change of small volumes of the fluid, 
given by 

9 = ^ = w(v + 2 f)' < 27 > 

where dot and prime denote the time and radial derivatives respectively. The 
components of field equations for spherically symmetric interior spacetime are 
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of the form 
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( 2 . 10 ) 

A(k 

W 2 

( 2 . 11 ) 


The dynamical equations are important in establishment of the instability 
range of collapsing stars. Misner and Sharp mass function furnishes the 
total amount of energy in a spherical star of radius ”Y” and facilitates in 
formulation of dynamical equations, given by jJ5] 


m(t, r) 


Y_( 

~2 y + W 2 ~^ 2 )' 


( 2 . 12 ) 


The matching conditions of adiabatic sphere on the boundary surface of 
exterior spacetime results from continuity of differential forms as m 

M = m{t,r ) (2.13) 


The dynamical analysis can be established by using conservation laws, we 
have taken conservation of Einstein tensor because the energy momentum 
tensor bear non-vanishing divergence in f(R, T) gravity. The contracted 
Bianchi identities imply dynamical equations that further leads to collapse 
equation, given by 


G u fV u = 0, 


;v 


G uv y 

^;v Ah 


= 0 , 


(2.14) 
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Bianchi identities in account with the Eg. (12.71) become 


P + P 


[1 + /t] © 



+ [1 + fr\ 


X Y 
pr ^ + 2 P±y 


+Zi(r,t) — 0, 

f W' (W Y' f' D \ 
(r + Pr)f^(l+fT){ti + P W +Pr( W + 2 Y -fy 

- 2 P±y } + h (p‘ - Y) + Z 2 (r, t ) = 0, 


(2.15) 

(2.16) 


where Zi(r,t) and Z 2 (r, t) are the corresponding terms including dark mat¬ 
ter components provided in Appendix as Eqs. (16 .1 H and (16.21) respectively. 
These equations are useful in the description of variation from equilibrium 
to the evolution. 


3 /(i?, T) Model and Perturbation Scheme 

The f(R, T) model we have considered for evolution analysis is 

f(R,T) = R + aR 2 + XT, (3.17) 

where a and A corresponds to the positive real values. Generally, a viable 
model represents the choice of parameters whose variation shall be in ac¬ 
cordance with the observational situations [[47]. Astrophysical models are 
selected by checking their cosmological viability, that must be fulfilled to ex¬ 
tract consistent matter domination phase, to assemble solar system tests and 
stable high-curvature configuration recovering the standard GR. The model 
under consideration is consistent with the stable stellar configuration because 
second order derivative with respect to R remains positive for assumed choice 
of parameters. 

The field equations in f(R, T) are highly complicated, their general so¬ 
lution is a heavier task and has not been accomplished yet. Evolution of 
linear perturbations can always be used to study the gravitational modifica¬ 
tions by avoiding such discrepancies. The concerned collapse equation can 
be furnished by application of linear perturbation on dynamical equations 
along with the static configuration of held equations leading to the insta¬ 
bility range. The dynamical analysis can be anticipated either by following 
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fixed or co-moving coordinates i.e., Eulerian or Lagrangian approach respec¬ 
tively El. Since universe is almost homogeneous at large scale structures 
that is why we have used co-moving coordinates. 

Initially all physical quantities are considered to be in static equilibrium so 
that with the passage of time these have both the time and radial dependence. 
Taking 0 < e -C 1 the perturbed form of quantities along with their initial 
form can be written as 


W{t,r) 

= 

Wo(r) + sD(t)w{r ), 

(3.18) 

X(t,r) 

= 

X 0 (r) + eD(t)x(r), 

(3.19) 

Y(t,r ) 

= 

Y 0 (r) + eD(t)y(r), 

(3.20) 

p(t,r) 

= 

Po(r) +ep(t,r), 

(3.21) 

Pr(t,r) 

= 

Pro(r ) + ep r (t,r), 

(3.22) 

P±(t,r) 

= 

P±o(r) + ep]_(t,r), 

(3.23) 

m{t , r) 

= 

m 0 (r ) + £fh(t, r), 

(3.24) 

R(t,r) 

= 

R 0 (r) + eD(t)ei{r)i 

(3.25) 

T(t,r) 

= 

T 0 (r) + eD(t)e 2 (r), 

(3.26) 

f(R,T ) 

= 

[Ro(r) + uRl(r) + AT 0 ] + eD(t)ei(r)[l 



+ 

2aR 0 (r)\ + r]D(t)e 2 (r), 

(3.27) 

In 

= 

1 + 2aR 0 (r ) + e2aD(t)ei(r), 

(3.28) 

fr 

= 

A, 

(3.29) 

0(t,r) 

= 

£0. 

(3.30) 


Without loss of generality, we have taken the Schwarzschild coordinate lo( r ) = 
r and apply perturbation scheme on dynamical equations i.e., Eqs. (12.15p and 
(12.161) , the perturbed Bianchi identities turn out to be 


P + 


2ep 0 


+ Xi 


1 T 2aR 0 

+ (1 + 2aRo)Zi p ] D — 0 , 

. f _, W' Q 
1 \ Pr + P W 0 +Pr 


2 y x 

— (Po + 2pj_o) + -rr(po + Pro ) 
r A 0 


K + l_ 

Wo r 


2aR' r 


0 


2pl| 


+2e—2 - — R! q ) + Xl{ 1 + 2aR 0 ) 
Aq a 0 


1 + 2aRo j r 

e 


+ 2a 


W 0 2 


(3.31) 

(e' 


Xq (1 + 2aR 0 ) 


D + D [Ai[(p 0 













. . w . . ,y. n - 2a 

+Prt0 W - 2 (p r „ +p±0 )(-) ] + A p- i - — 


+PrO 


^0 


2aR' n 


r Wo 1 + 2aRo 
+ (1 + 2aR 0 )Z 2p \ = 0. 


)) } + ^ + e ^° ~~ 


Al [ Pr0 + po Wo 

2 aR' n 


1 -)- 2aRo 


(3.32) 


For the sake of simplicity we assume that e\ = e 2 = e and set Ax = A + 1, 
Zip and Z 2p are provided in appendix. 

Elimination of p from Eq. (I3.3ip and integration of resultant with respect 
to time yields an expression for p of following form 


P = ~ 


2ep 0 


1 + 2aRo 


2 y 


x 


+ Ai — (po + 2pj_o) + ~^(Po + Pro) f + (1 + 2aRo)Zi p 


D. 


(3.33) 

The perturbed on field equation Eq. (l2.1ip leads to the expression for p±, that 
turns out to be 


P± = 


D f (1 + 2aR 0 )y 


VF 0 2 


D, 


- 2ae ^ + 


A 


2ae 


P±0 \i P °Jl + 2aRo 


(3.34) 


effective part of the field equation is denoted by Z 3 , given in appendix 
Eq. (16.51) . Matching conditions at boundary surface reveals 


^ a ^ n 

Pr = 0, p± = 0 


(3.35) 


Above equation together with perturbed form of Eq. (12.111) can be written in 
the following form 

D(t) — Z±(r)D(t) = 0, (3.36) 

provided that 


= 


tWq 


(1 + 2aRo)y — 2 aer [1 + 2aRo 


2ae 


P±o + A <! — (po + 2p ±0 ) 


x 


+ — {po + Pro ) f + (1 + 2ai? 0 )Zi p + — 

-&o J M 

The solution of Eq. fl3.36p takes form 

D(t) = -e^. 


(3.37) 


(3.38) 
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The terms appearing in Z 4 are presumed in a way that all terms remain 
positive to have a valid solution for D. Expansion-free condition and stability 
range is discussed in the following section. 


4 Expansion-free Condition with Newtonian 
and Post-Newtonian Limits 


The models with an additional zero expansion condition delimitate two hy¬ 
persurfaces, one separates the external Schwarzschild solution from the fluid 
distribution and other is the boundary between fluid distribution and internal 
cavity. Such models have extensive astrophysical applications where cavity 
within fluid distribution exists and are significant in investigation of voids at 
cosmological scales [42]. The spongelike structures are termed as voids ex¬ 
isting in different sizes i.e., mini-voids to super-voids mm accompanying 
almost 50% of the universe, considered as vacuum spherical cavities within 
fluid distribution. 

Implementation of linear perturbation on Eq. (l2.7p and Eq. fl2.12p respec¬ 
tively, implies 


e = i(^ +2 ? 


hi 


VEn V X, 


m 0 = - 1 - 




m — — 


D 

Xl L 


r l y - 


x 

Xn 


2^y 


+ (i - xt) 


(4.39) 

(4.40) 

(4.41) 


The expansion free condition implies vanishing expansion scalar i.e., 0 = 0 
implying 


x 

Xn 


= -2 


y 


(4.42) 


In order to present the dynamical analysis in Newtonian (N) and post- 
Newtonian (pN) limits, we assume 


Po ^ Pro , Po P±o- 


(4.43) 


The metric coefficients upto the pN approximation in e.g.s. units are taken 
as 


W 0 = 


1 - 


Gmo 

9 ’ 

c z r 


X 0 = 


1 + 


Gm 0 

9 ’ 

c z r 


(4.44) 
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where c denote speed of light and G stands for gravitational constant. Ex- 
pression for ^ can be obtained from Eq. (14.401) as 


K = ~ m o 

X 0 r(r — 2m 0 ) ’ 

Eq. il 1. Kill together with (12. 10p implies 


(4.45) 


K 

hho 


—---- v [^(AlPrO + Vo - Ro 

2 r(r — 2m 0 )(l + 2 an 0 + ran 0 ) L 

— 3a/2g) + 2ar(R 0 — 2rR' 0 + 4:R' 0 m 0 ) + 2m 0 ] . (4.46) 


The static configuration of first Bianchi identity is identically satisfied while 
second provides a fruitful result in terms of dynamical equation. Substitution 
of Eqs. (14. 45 ft and (j4.46[) in statically configured Eq. (j2. 16[) and after some 
manipulation the dynamical equation in relativistic units yield 
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+2ar(i?o — 2r.Ro + 4i?omo) + 2mo 


— [r 3 (Aip r0 + Apo — #o 


Po + PrO + 


2ai?o(r 


3ai?o) 


2m n 


x 


/ 3m 0 1 r 3/\ 
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3ai?o) + 2ar(i? 0 — 2ri?' 0 + 4/?o?Ti 0 ) + 2m 0 — Rq] }]] • (4.47) 
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In c.g.s. units, we may write above equation as 


PrO 
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2 
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+ W 


(4.48) 


The terms of order c° and c -2 belongs to N and pN-approximation respec¬ 
tively. One can expand Eq. (14.48j) upto c -2 and separate the terms of N and 
pN limits to distinguish the physical quantities lying in various regimes. 

The use of expansion-free condition in Eq. (I3.33P modifies p to following 
form 


P = ~ 


-r~ o“"TT + Ai— ((Pro ~ P± o)) + (1 + 2otRo)Zi p 
1 + 2aR 0 r 


D. (4.49) 


The Harrison-Wheeler type equation of state describing second law of ther¬ 
modynamics relates p and p r in terms of adiabatic index T as 

p r = T Pr0 p. (4.50) 

PO + PrO 

T measures the fluid’s compressibility belonging to its stiffness. Inserting p 
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from Eg. (14.4911 in Eq. fl4.50p . we have 


Pr 


p PrO 
P0 + PrO 


2ep 0 
1 -|- 


. \ 2 y 
+ X - {{PrO 


P± o)) + (1 + 2ai?o)^ip 


D. 


(4.51) 

In view of dimensional analysis, it is found that terms of p r and Poy^ lies 
in post post Newtonian (ppN) era and thus can be ignored from the terms of 
N and pN approximation. Since we are going to exclude p r , so it is intuitively 
clear that instability range is independent of T, no compression is introduced. 
Use of expansion-free condition together with the expression found for D, it 
follows that 


2 aZ 4 


1 f e'+ 2e^ + 2^-R' 0 ) + X 2 (l + 2aR 0 ) 
A 0 r 




Xq (1 + 2ct Rq) 


2 a 


+Allp " ( iy ) ' " 2(Pr0 + P±o) X ] + V ~ 1 + 2aR„ 


X PrO + PrO 


2aK 


1 + 2 aR 


;)) } + A (e' + e[g - j + l 1 + 2aI >o) z 2 r ~ ». 

(4.52) 


For simplification of above expression, we take relativistic units and assume 
that po Pro, Po P±o- O n substitution of expressions for Z 4 ,Wo,X 0 , 
and Z 2 p respectively from Eqs. (I3.37p . (I4.44p . (14.45(1 and ((6.4(1 yields a very 
lengthy expression defining the factors affecting the instability range at N 
and pN limits. The expanded version of Eq. (14 .52 (1 is large enough, therefore 
we are quoting only the results obtained from the collapse equation together 
with the restrictions to be imposed on physical parameters. It is clear from 
Eq. (14.47(1 that p' r0 < 0, provided that all the terms maintain positivity to 
fulfill stability criterion. Negative values of p' r0 depicts decrease in pressure 
with time transition, leading to collapse of gravitating star. Furthermore, 
using the e.g.s. units, it is found that the terms of p r and Pof^ does not 
take part in evolution for N and pN approximations since these terms belong 
to ppN limit. The analysis of terms lying in N and pN limits imply few 
restrictions to be imposed on physical quantities for discussion of instability 
range. These are listed as follows: 


• Newtonian Regime: The constraints on material parameters are 


Pro > P±o, a 2 rR' 0 < 1 + 2 aR 0 , 


2 aR' 0 

1 T 2cti?o 


> Po - e ' ■ 
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The gravitating body remain unstable as long as the inequalities hold 
in N approximation. 


• post-Newtonian Regime: In pN limits following restrictions are 
found to execute the instability range 



r 


Pro > Rlo, r > 2m 0 , 


r + mg 



5 Summary and Results 

The mysterious content named as dark energy (DE) occupying the major 
part of universe is significant in the description of cosmic speed-up. The 
modified gravity theories are assumed to be effective in understanding cos¬ 
mic acceleration by induction of so-called dark matter components in the 
form of higher order curvature invariants. Among such theories, f(R, T) 
represents non-minimal coupling of matter and geometry. It provides an al¬ 
ternate to incorporate dark energy components and cosmic acceleration {52j . 
Thus consideration of f(R, T) for dynamical analysis is worthwhile, covering 
the impact of higher order curvature terms and trace of energy momentum 
tensor T. This manuscript is based on the role of viable f(R,T) model in 
establishment of instability range of spherically symmetric star. 

Our exploration regarding viability of the f(R, T) model reveals that the 
selection of f(R , T) model for dynamical analysis is constrained to the form 
f(R, T ) = f(R) + AT, where A is arbitrary positive constant. The restriction 
on f(R, T) form originates from the complexities of non-linear terms of trace 
in analytical formulation of field equations. The f(R , T) form we have chosen 
mainly is f(R, T) = R + aR 2 + AT, in agreement with the stable stellar 
configuration and satisfies the cosmic viability. The matter configuration 
is assumed to have unequal stresses i.e., anisotropic with central vacuum 
cavity evolving under expansion-free condition. Zero expansion condition on 
anisotropic background reveals the significance of energy density profile and 
pressure inhomogeneity in structure formation and evolution. 

The field equations framed in f(R,T ) gravity are formulated and their 
conservation is considered to study the evolution. Conservation laws yield 
dynamical equations that are significant in formation of collapse equation. 
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In order to examine the variation from static equilibrium, we introduced 
linear perturbation to all physical parameters. The expressions for perturbed 
configuration of field equations reveal expressions for energy density p and 
tangential pressure p±. The second law of thermodynamics relating radial 
pressure and density with the help of adiabatic index T is considered to 
extract p r . 

On account of zero expansion it is found that fluid evolution is inde¬ 
pendent of T, rather instability range depends on higher order curvature 
corrections and static pressure anisotropy. Recently, the dynamical analysis 
of isotropic and anisotropic spherical stars in f(R, T) has been studied in 
[3S» HQ]. It is found that perturbed form of dark source terms of collapse 
equation also has the contribution of trace T, affecting the stability range. 
Thus non-minimal coupling of the higher order curvature terms and trace of 
energy momentum tensor imply a wider range of stability, however, the fluid 
evolving with zero expansion might cause drastic and unexpected variations. 
As expansion-free condition produces shear blow-up in gravitating system, 
so it is very captivating to extend this work for shearing expansion free case. 
The results are in accordance with [13] for vanishing A, for vanishing a and 
A corrections to GR solution can be found. 

In addition to model fj3. 17j) . nature of various f(R, T) models i.e., f(R, T) = 
R + aR n + AT, f(R, T) = R + aR 2 + £ + XT and f(R, T) = R+^ +XT 
has been briefly discussed in this section, as follows: 


• f(R, T) = R + aR n + AT: The model f(R, T) = R + aR n + AT is 
viable for any n > 2 and positive constants a and A. The collapse 
equation for such model with zero expansion is of the form 


W' 

A p + Ai 'x pr + (p + Pr)j7T + Pr ( — 


Wn 


an(n - l)Ro R ' 0 


1 + anR, 


>n—1 


2 _ 

-P± 

r 


+D 

+ ^3 p, 


w' 2 . . an(n - 1 )(Rq 2 e)' 

AlfPo + Pro) 777- + Al — \PrO ~ P± o) “ (A + AiP r0 )--- 

Wq r 1 + anRp 


(5.53) 


where pressure stresses p r ,p± can be generated from perturbed held 
equations and perturbed energy density is 


P 


anin 


1)(R { 


n —2 . 


1 + anR 0 


g) oc 

Po + Ai(p r o — p±o)~rr + Z‘ 


in—1 


An 


4 p 


D, 


(5.54) 
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where Z 3p and Z 4p depict the perturbed dark source terms. The N and 
pN limit of this model reveals that the term p r and belong to 

ppN limit and so do not contribute in evolution. In Newtonian limit 
the physical quantities must satisfy the following conditions 


Pro > V -an a 2 rn(n — l)i?Q 2 R' 0 < 1 + naR g h 
nan(n — 1 )Rq~ 2 R' 0 


-)Tl— 1 


> Po ~ 


1 + HOiRq 

The constraints on physical quantities in pN regime are 

r > 2mo, - (xn(n — 1) RT~ 1 RL H--) < e' + A, 

r + mo r 


2ae-(l + na R S ~^ > 


er 


(r + m Q ) J 


(r — 2m 0 ) > 


(n - 1 )KpK 

™-l 

K 0 


,Pr 0 > RLO- 


The dynamical analysis of various models involving higher order cur¬ 
vature terms, combined with the trace of energy momentum tensor can 
be presented for n > 2. 

• f(R,T) = R + aR 2 + ^ + XT and f(R,T) = R + ^ + XT: The 
perturbed form of Bianchi identity for f(R, T) = R + aR 2 + % + XT, 
where /i is arbitrary constant leads to the expression for p as follows 


P = 


2ep 0 


x 


^~ 2 Al 1 j^ _ (P r0 — P-Lo) \ + (1 + 2qRq — fi R 0 )Z\ p 


(5.55) 


[ 1 + 2aRo — p 4 Ro 
The evolution equation becomes 
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+Pr0 {r + W 0 l + 2aRo-^Ro 2 
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+ A ^e' + e[po 
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, ( , W' 

Al [ Pr0 + p 0 W 0 

2 aR' 


1 + 2aR 0 


(5.56) 


16 



Z$ p denote the perturbed dark source entries. The N and pN limits are 
obtained by avoiding the terms lying in ppN region. To maintain the 
viability of model in Newtonian era, following inequalities must hold 


c,R' a +^RfK„ < l+2 a Wfl„- 2 , + o > 

1 + 2 aRo — fJ- R 0 

In pN regime the system remains stable as long as following ordering 
relations are satisfied. 

- (cti?g + 3 R' 0 ) < e' + AT 0 , 

r + m 0 

2ae - ^Rf - (1 + 2 aR 0 - ^Rf)^ > ^ ~ ^ 

4 

The collapse equation for f(R, T ) = R + ^ + AT can be obtained 
by setting a = 0 in Eq. (I5.56I) , likewise the restrictions on physical 
quantities can be found. 
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6 Appendix 
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